function moment = hwozd_mi(I,p,q)
    xs = (I.x - I.x0) .^p;
    ys = (I.y - I.y0) .^q;
    moment = sum(sum( (xs .* ys) .* I.image));
